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Abstract 

We study the translocation of a flexible polymer through a pore in a 
membrane induced by its adsorption on trans side of the membrane. When 
temperature T is higher than T c , the adsorption-desorption transition temper- 
ature, attractive interaction between polymer and membrane plays little role 
in affecting polymer conformation, leading to translocation time that scales 
as r ~ 1? where L is the polymer contour length. When T < T c , however, the 
translocation time undergoes a sharp crossover to r ~ I? for sufficiently long 
polymers, following the second order conformational (adsorption) transition. 
The translocation time is found to exhibit the crossover around T = T' c , which 
is lower than T c for polymers shorter than a critical length(iV < N c ). 
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I. INTRODUCTION 



Biopolymers exhibit various transmembrane conformational states determined by their 
intrinsic properties as well as by the interactions with the surrounding environments. For 
example, polymer flexibility, capability of hydrogen bond formation between segments, and 
hydrophobic interaction with membrane play crucial roles in conformation of membrane- 
bound proteins. An important aspect of membrane-protein interaction can be found in the 
dynamics of proteins. As the experimental observations of protein translocation are available 
only in limited situations]!! theoretical investigations are also in a primitive stage. However, 
a few studies show that polymer translocation across the membrane is significantly affected 
by the surrounding environments and membrane-protein interactions.!!! Theoretical studies 
have been focused on the driving mechanisms of polymer translocation across the membrane 
such as chemical potential asymmetry! and its random modulation^ chaperon bindings,!! 
and lipid molecular density imbalance.! It was shown that a minute chemical potential bias 
per segment of the chain can bring about a significant effect on translocation, due to the 
long-chain connectivity.!' Bi In this paper, we investigate the effect of attractive membrane- 
polymer interaction yielding adsorption transition as another possible driving mechanism. 

To highlight the effects arising from long chain nature, we consider a rigid plane of 
negligible thickness with a hole small enough to allow only a single segment passage, as 
a minimal model of membrane with a translocation pore as shown Fig. 1. In this model, 
polymer segments interact with the membrane attractively in the trans side, while those in 
the cis side recognizes the membrane only as an impenetrable wall. No specific interaction 
between the segments and pore is assumed. Furthermore we neglect the mutual interaction 



2 



between the segments because we are interested in the most dominant effects resulting from 
the chain connectivity. The attraction, if it is strong enough to bring about desorption- 
adsorption transition, would drive the polymer to the trans side as shown in Fig. 1 to 
minimize the free energy. 

The total free energy of the transmembrane polymer has two distinct contributions. One 
is the entropic contribution affected by the steric interaction with the membrane, and the 
other is the adsorption energy due to attraction in trans side. At high temperatures, thermal 
fluctuation dominates and the polymer conformation is nearly the same as the case without 
attractive interaction. As the polymer is situated in nearly symmetric environments in this 
case, it is driven to either side of the membrane with equal probability. At temperatures 
below the adsorption-desorption transition temperature, on the other hand, polymer con- 
formation in the trans side becomes collapsed(adsorbed) and the polymer is driven sponta- 
neously to the trans side to minimize the total free energy. Near the transition temperature, 
however, competition between thermal fluctuation and membrane attraction determines the 
chain conformation, while the translocation dynamics is governed by the interplay between 
thermal diffusion and segmental energy bias due to interaction with membrane. 

In section II, conformation of transmembrane polymer and the associated free energy are 
derived as a function of temperature and interaction potential parameters. In section III, 
we consider the dynamics of polymer translocation as a stochastic process crossing the free 
energy barrier, where the effects of attraction and temperature on translocation dynamics 
are discussed. Section IV gives summary and conclusion. 
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II. STATISTICAL MECHANICS OF A POLYMER ONE-END- ANCHORED ON 

ATTRACTIVE SURFACE 



The transmembrane polymer can be regarded as two independent chains anchored at the 
spot of pore(Fig. 1). In this section, we consider the chain conformation and free energy 
of one-end-anchored chain to calculate the free energy of a transmembrane polymer. The 
Green's function of a polymer G(r, r'; n) is defined as the statistical weight of an n-segment 
polymer whose end positions are located at r and r' respectively. For a polymer with Kuhn 
length b, G(r, r'; n) satisfies the following differential equation: 
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G(r, r'; n) = 8(r — r'). 
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In unbound free space it is reduced to a familiar Gaussian form: 
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For a polymer subjected to an attractive wall(z = surface), the attractive potential 
effect can be incorporated in terms of the following boundary condition,@ 
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This boundary condition consideration is valid only for the case that the attractive inter- 
action range is short enough compared with other length scales involved such as polymer 
radius of gyration or adsorption layer thickness. The c is a convenient parameter descriptive 
of surface interaction the details of which this mesoscopic description usually bypasses. Al- 
though, for flexible chain, the Edwards equation is more adequate to describe the adsorbed 



state incorporating details of surface interaction, this mesoscopic description, also applicable 



to semiflexible chains, can easily yield useful results for a full range of temperatures. For 
the cases of weak adsorption it <C ksT and flexible chain a>ft considered throughout this 
paper, we will use the Edwards equation to determine the interaction parameter c in terms 
of temperature, attraction strength(wo) and range(a). 

The Green's function in half space(z > 0) can be written as 



G(r, r'; n) = G (x, x'; n)G (y, y'\ n)G(z, z'\ n) 



(4) 



where Go(x, x';n) and Go(y,y';n) are one dimensional Gaussian propagators analogous to 
Eq. |2|, and the z component is given bylE 
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where erfc(w) = (2/ \/7r) exp(—u' 2 )du' is the complementary error function, and 

w = (n/6) 1/2 bc. 

The z component of Green's function has the following limiting behaviors: 
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as it should be. The latter is the ideal chain Green's function in the presence of absorbing 
boundary at z — 0. Now that \w\ is defined as the ratio of free polymer size(~ n 1 ^ 2 ) and 
adsorbed layer thickness(~ |c| _1 ) for low temperature adsorbed states(w < 0, c < 0), it 
can be regarded as a parameter describing the effects of surface attraction or temperature. 
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For high temperature desorbed states(u> > 0), w measures the strength of steric repulsion, 
which also entails the surface attraction effect. 

The conformational partition function of a polymer with one end anchored on attractive 
surface is given by0 
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z'-*0 Jo 



exp(iy 2 )erfc(w) 
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1 - 2n- 1 ^w, w ~ 
2exp(u> 2 ), w —1, 

where the bulk value of partition function q n is omitted throughout this paper, where q is 
the partition function of a single segment. To see the conformational states of anchored 
polymer explicitly, let us see the conditional probability distribution of free end defined by 



Pn(z\z' = 0) 



Q(n 



-G{z, z' = 0;n), 



(12) 



which has the following limiting forms: 

{nb 2 /3)- 1 zexp(-3z 2 /2nb 2 ), w > 1 
P n (z\z' = 0) = \ (Trnb 2 /Q)~ 1 / 2 exp(-3z 2 /2nb 2 ), w = (13) 
|c| exp(— \c\z), w <C — 1. 

For w 3> 1, P n (z\z' = 0) is equal to the distribution without attraction, which corresponds 
to the high temperature desorbed phase. On the other hand, for w —1, P n (z\z' = 0) 
represents the distribution for low temperature adsorbed phase of a polymer subject to a 
surface of sufficiently short range attraction. The distribution at the adsorption-desorption 



transition point corresponds to the case w = 0. 



Using the partition function of a one-end-anchored polymer in Eq. [10], its free energy is 
given by 



F(n) = -k B T\ogQ{n) 



-k B T 



w 2 + logerfc(w) 



(k B T/2) log n + const, w ^> 1 
2k B Tn^ 2 bc(6n)-^ 2 , w~0 
—kBTnb 2 c 2 /6 + const, w -C — 1. 



(14) 
(15) 



(16) 



Note that the F(n) changes sign at w — 0, which corresponds to adsorption-desorption 
transition point. At high temperature (w 3> 1), we can see the polymer free energy, the 
value relative to that of free polymer, goes like logn due to the anchorage effect. At low 
temperatures(u> —1); adsorption free energy scales as n involving segmental free energy 
as -k B Tb 2 c 2 /Q. 

Since the transmembrane polymer is composed of two independent one-end-anchored 
chains with length n and iV — n in cis and trans side respectively(Fig. 1), its partition 
function is given by 



Q(n) = Q cis (N - n)Q trans (n, c) 



(17) 



where the effect of attraction appears only in the trans side partition function in terms of c. 
This partition function has the following limiting behaviors: 
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The free energy as a function of translocation coordinate n, defined as the number of segments 
in trans side, is, apart from additive constants, 



T(n) = — k B T log Q(n 
1 
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T > T c 

\k B T\og{N -n), T = T C (21) 

\k B T\og{N - n) - k B T^f, T < T c . 
In Fig. 2, this free energy function is depicted for different values of c. T{n) exhibits nearly 
symmetric barrier for positive c(T > T c ), while it does exhibits a linearly slant shape for 
negative c(T < T c ). This signifies that the polymer is driven to the trans side due to the ad- 
sorption for c < 0. One remarkable point here is that J-{n) slants to the trans side at T = T c 
even though the polymer in the trans side is not in its adsorbed state. The n segments in 
the trans side feel neither steric interaction nor attraction. The cancelation of these two 
effects makes the total free energy be determined only in terms of the N — n segments in 



the cis side. 



Determination of c 

In the previous section, the free energy of a transmembrane polymer is determined in 
terms of interaction parameter c, which is defined as the logarithmic derivative of polymer 
Green's function at the surface (Eq. |^). To determine c in its adsorbed phase in terms of 
interaction potential parameters, we consider a microscopic model of flexible polymer weak 
adsorption in this subsection. As the polymer conformation in its adsorbed phase is not 
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affected by the anchorage, we will consider the adsorbed polymer without anchorage for 
simplicity 

Let us choose the interaction potential between a monomer and the membrane as 

—Mo, < z < a, 

(22) 

0, z>0, 

where uq > is the strength and a is the range of attraction. Under the conditions uq <C ksT 
and a ^> 6,0 the polymer Green's function obey the Edwards equation 

d 

G{z, z'; n) = ~^~G(z, z 1 ; n). (23) 
It can be expanded as 

G(z, z'; n) = £ exp^nA^^)^') (24) 

i 

where Aj and ipi(z) are the eigenvalues and eigenf unctions of the Schrodinger-like equation, 

^i{z) = \^i{z), (25) 
subject to the boundary condition 

il>i{z = 0) = 0. (26) 

By solving this eigenvalue problem, we can get the adsorption-desorption transition temper- 
ature T c and the value of c as a function of T, a and u for adsorbed state. 

For T < T c , the ground state dominance approximation, G(z, z';n) ~ 
exp(— n\ )tpQ(z)tpo(z') , works well for a long chain, where A is the lowest eigenvalue and 
ipo{z) is the corresponding eigenfunction of Eq. The ground state eigenfunction in z 
direction is given by 
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where k = {Q/b 2 ) 1 ' 2 ^ + A ) 1/2 and re = (6/6 2 ) x / 2 1 A 1 1/2 with A < for T < T c . The 
transition temperature T c can be determined as 



k B T c 
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below which a bound state of negative eigenvalue exists. The continuity condition of ipo(z) 
and ip' (z) at z = a gives the following relation: 



Ka = —ka cot(ka) 



(29) 



which determines Aq for T < T c . The interaction parameter c for T <T C can be defined as 



c = 



1 d 

i[)q(z) dz 
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This microscopic definition of c can be identified with the previous definition if the interaction 



is short-ranged. From the solution of Eq. 29 and Eq. 30, c can be written as the following 



functional form for T < T r : 



c = —a 1 H(a) 



(31) 



where a = (a/&)(6/?wo) 1//2 , and H(a) is depicted in Fig. 3. As T approaches T c from below, 
a goes to 7r/2 and H(a) goes to zero, which yields 



8a T r 
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for T < T c . At sufficiently low temperatures(T <C T c ), on the other hand, a ^> 1 and 
H(a) ~ a. In this regime, c can be approximated as 
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c~ -b-^Puo) 1 ' 2 ~ -T- 1 / 2 , (33) 

which is independent of interaction range a. Eq. ^ and Eq. |33| are useful to examine the tem- 
perature dependence of polymer descriptions derived from the boundary condition(Eq. 
which will be demonstrated in the next section. 

III. DYNAMICS OF POLYMER TRANSLOCATION 

As shown by us,! the translocation of a polymer can be thought of as a one dimensional 
diffusive process of the translocation coordinate n crossing the effective potential !F(n). The 
probability density of n(t), P(n,t), is described by Fokker-Planck equation 

^P{n,t) = C FP {n)P{n,t), (34) 

where C F p(n) is the Fokker-Planck operator given by 

C FP {n) = __D(n)exp[-/3^(n)] — exp[/3^(n)], (35) 
tr on on 

with D(n) defined as the whole chain diffusion coefficient. The translocation time of a 
polymer can be defined as the mean first passage time r(n; no), the time for diffusion from 
n to n, which satisfiesEl 

C FP (n )r(n;n ) = -1, (36) 

where £ FP (no) is the backward Fokker-Planck operator given by 

Id d 
C FP (n ) = -n: explpFtno)}—- r>(n ) exp[-/3.F(no)]^— • (37) 
tr on On 

Using the following boundary conditions 
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-r(n;n = 1) = (38) 



dn 

r(n;n = N- 1) = 0, (39) 

the solution of the above backward equation(Eq. |36|) can be obtained and the polymer 
translocation time is determined as 



T = T 



(N - 1; 1) (40) 



f N 1 dn 1 e PHn) r dn ' e -W). (4i) 
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By putting the diffusion coefficient D(n) to be a constant Z) for simplicity, the translocation 
time as a function of w is shown in Fig. 4, where a drastic change in r occurs near w = 0. 
Because w or c is expected to be a monotonic function of temperature as shown in Eq. [IT] 
and Fig. 3, Fig. 4 can be understood as a sharp dynamic crossover induced by temperature 
change. Note that r ~ \w\~ 2 ~ |c|~ 2 for T T c , while r changes smoothly only by a small 
amount for T > T c . 

The Fig. 5 is a plot of r, which is obtained using the free energy(Eq. ^0|) with c given by 
Eq. |3l], as a function of iV for four different temperatures ranging from T ^> T c to T T c . 
For T ^> T c , i.e., /3uq ~ 0, the transmembrane polymer does not feel attractive interaction 
and translocation time(Eq. (41]) is reduced to 
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2D 

where L = Nb is the contour length of the polymer and the diffusion coefficient is given 
by D = ksT/lN'y), with segmental hydrodynamic friction coefficient 7 following the Rouse 
model.0 In this high temperature regime, translocation time is proportional to L 3 (Fig. 5, 
A), which is the same result as the case without the attractive interaction.! 
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As the temperature is lowered, the free energy curve tends to slant down to the right 
and the prefactor of r decreases smoothly until T = T c as shown in Fig. 4. At T = T c , the 



translocation time becomes 



2D 



(43) 



where the prefactor is reduced from 7r 2 /8 to 2/3. In this case, the steric interaction in 
the trans side is canceled off by the attractive interaction, which drives the polymer to the 
trans side, but r is still proportional to L 3 (Fig. 5, B). Assuming D increases as T increases 
like D ~T a with a > 0, 



r(T > T c ) 
r(T = T c ) 



3tt 2 (T c , 
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Thus, as can be seen in A and B of Fig. 5, the r increases as T decreases to T c , which results 
from the suppression of thermal diffusion. The r attains its maximum at T = T c , since for 
T < T c , the polymer adsorption speeds the translocation as discussed below. 

For T <^ T c , c behaves as c ~ — T _1 / 2 (Eq. |S3|) and the translocation time(Eq. is 
given by 
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which is depicted in Fig. 5, D. In this low temperature regime, the free energy of transmem- 
brane polymer is dominated by the attraction energy which is linearly proportional to n, 
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where the prefactor of n can be interpreted as the energy per segment in trans side, or trans- 
membrane segmental energy bias that features as "chemical potential difference" between 
the two sides of membrane.! 

We find that the crossover of scaling behavior from r ~ L 3 to r ~ L 2 sharply occurs 
around T = T c , which can also be inferred from Fig. 4. The conditions for the crossover 
can be analyzed in detail as follows using the Edwards equation approach for T < T c . As 
shown by us,! segmental energy bias approximately larger than ksT/N can bring about the 
crossover in translocation from r ~ L 3 to r ~ L 2 due to chain connectivity(cooperativity). 
For the adsorption- driven translocation when T < T c , the condition for the crossover can be 
written as 

|Ao| > p (49) 

where Ao, the lowest eigenvalue in Eq. |25], is the segmental free energy in units of ksT in 
adsorbed state. This condition(Eq. EST) for T near T c can be converted to 



T < T' C (N) = T c 



8a / 6 X 1 / 2 



(50) 



7T 2 \Nb 2 , 

where the relation Ao = — b 2 H 2 (a)/6a 2 (see Eq. |3l]) was used. This condition implies that the 
crossover occurs at T' c smaller than T c for finite chain length N. That is, for T' c < T < T c , 
t ~ L 3 although the transmembrane polymer is adsorbed in trans side. For a numerical 
example for the difference between T c and T c ', we consider the case, a/b = 10 and N = 60000, 
which yields T' c ~ 0.9T C . The C of Fig. 5 depicts this situation. For temperatures less than 
T c , the condition of Eq. H9I yields the critical chain length 



NAT) , ^ ) (81) 
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around which the crossover from r ~ L 3 to r ~ L 2 occurs. Indeed this N c is consistent with 
the crossover point in C of Fig. 5. 

Table 1 summarizes the conditions for the crossover and associated chain conformations 
for diverse temperatures and chain length. The meaning of two critical temperatures T c , T' c 
as well as the critical chain length N c for chain cooperativity, can be understood as follows. 
The segmental free energy Ao(T) vanishes at T = T c , while N\\q(T) \ becomes larger than 1 
for T < T' c or N > N c , which are the conditions for a cooperative, facilitated translocation 
driven by an energetic bias. In this study, we do not incorporate the solvent effect which 
will induce interaction between segments and thus affect the polymer free energy function 
significantly. Yet the crossover in r is seen to occur irrespective of excluded volume effect, 
as a consequence of the long-chain nature of the polymer interacting with the membrane.! 
Incorporation of the effect will modify only T c and T' c . 

IV. SUMMARY AND CONCLUSION 

We considered the dynamics of flexible polymer translocation as a diffusive process cross- 
ing the free energy barrier which we derived in section II. The translocation dynamics is 
found to be strongly correlated with transmembrane chain conformation, both of which are 
varied depending upon the competition between thermal fluctuation and attractive polymer 
interaction with the membrane. The polymer adsorption-desorption transition point (T c ) 
divides the length scaling regimes of translocation time, r ~ L 3 and r ~ L 2 , for sufficiently 
long chains(Y > N c ). However, for shorter chains(Y < N c ), a crossover of two regimes 
is found to occur at a lower temperature, T' c , although the conformational transition still 
occurs at the same transition temperature, T c . We used the mesoscopic theory of polymer 
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adsorption using the machinery of boundary condition, supplemented by the microscopic 
Edwards equation approach to consider the case of flexible chain weak adsorption in detail. 
This theory can be adapted to the more biologically-relevant situations of semiflexible chain 
strong adsorption, if the interaction parameter can be appropriately given. 
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FIGURES 

FIG. 1. Translocation of a polymer through a pore. Polymer segments are adsorbed in the 
trans side due to the attractive interaction at sufficiently low temperatures. 

FIG. 2. Free energy T{n) as a function of translocation coordinate n. (A. be = 1.0, B. be = 0, 
C. be = -0.5, D. be = -1.0) 

FIG. 3. H(a) versus a = (a/&)(6/3« ) 1/2 - At T = T c (a = tt/2), H(a) = and 
c = -a" 1 F(a) = 0. For T <C T c (a > 1), fl"(a) ~ a and c ~ -6~ 1 (6/3u ) 1/2 (dotted line). 

FIG. 4. Translocation time r normalized by ro = L 2 /(2D) versus \w\, with u; = bciN/6) 1 / 2 , 
(A. c> 0(T > T c ), B. c < 0(T < T c )). Note that r/r ~ w~ 2 for T <C T c . Inset: The same figure 
in linear scale, where sharp transitional behavior occurs at c = 0(T = T c ). 

FIG. 5. Translocation time r in units of n = ~/b 2 /uo versus chain length AT, where 
a = 106 < Nb, (A. T = 10T C , B. T = T c , C. T = 0.9T C , D. T <C T c ). r ~ L 3 for A and B, 
and r ~ L 2 for D. Crossover between them can be seen for C near N = N c as the chain length 
varies. 
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Table 1. Length scaling behavior of polymer translocation time r and associated chain 
formations as a function of temperature and chain length. 





shorter chain(iV < N c ) 




longer chain(iV > A^ c ) 


T>T C 


t ~ L 3 (desorbed) 


T>T C 


t ~ L 3 (desorbed) 


T' C <T < T c 


r ~ L 3 (adsorbed) 


T <T C 


r ~ L 2 (adsorbed) 


T < T' c 


t ~ L 2 (adsorbed) 
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Fig. 1, Park, Polymer Translocation 



Fig. 2, Park, Polymer Translocation 




Fig. 3, Park, Polymer Translocation 




Fig. 4, Park, Polymer Translocation 




Fig. 5, Park, Polymer Translocation 




